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FOREWORD 


This report, one of a series of four, describes ana- 
lytical and experimental results for flow in a convergent- 
divergent nozzle with massive wall blowing. The work was 
supported by the Nuclear Systems Division, NASA Lewis 
Research Center, under Grant NGR 03=002-213. Mr. Albert 
F. Kascak was the technical manager. 
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SUMMARY 


An analytical and experimental investigation has been 
conducted to determine the effect of massive wall injection 
on the flow characteristics in a nozzle. The experiments 
were performed on a water table with a porous-nozzle test 
section. This had 45° and 15° half angles of convergence 
and divergence, respectively, throat radius of 2.5 inches, 
and throat width of 3 inches. The hydraulic analogy was 
employed to qualitatively extend the results to a compres- 
sible gas flow through the nozzle. 

An analysis of the water-table flow was made using a 
one-dimensional flow assumption in the continuity and 
momentum equations. Allowance was made for fluid injec- 
tion along the nozzle wall. Predictions of the flow 
conditions in the nozzle were made for injection rates up 
to 30% of the inlet flow rate. 

Experimental results from the water table include 
contours of constant Froude and Mach number with and with- 
out wall injection. Photographs of the visualize! flow 
showing the injected layer are also presented for a range 
of injection rates. 

An analysis of a compressible flow in a nozzle was 
made in a manner analogous to that for the water flow. 

It is shown that the effect of blowing is to move the 
sonic position downstream of the geometric throat. Sim- 
ilar results were determined for the incompressible water- 
table flow. 

Limited photographic results are also presented for 
the injection of air, CO 2 > and Freon-12 into a main-stream 
air flow in a convergent-divergent nozzle. Schlieren 
photographs were used to visualize the flow, and qualita- 
tive agreement between the results from the gas tunnel 
and the water table is good. 

All the results are discussed in detail. Conclusions 
are presented which should be an aid in the preliminary 
design of a rocket nozzle for use in the gas core nuclear 
rocket . 


INTRODUCTION 


Advanced rocket engines, in particular the gas core 
nuclear rocket, are expected to produce gas temperatures 
at the entrance to the nozzle in the neighborhood of 
20,000° K (1). Such temperatures are orders of magnitude 
above those which can be tolerated by the structural com- 
ponents of the nozzle fabricated from existing materials. 
Therefore, the nozzle wall must be thermally protected to 
a much higher degree than that required by the current 
generation of rocket nozzles. 

A major portion of the nozzle heating is due to the 
thermal radiation from the high temperature gases in the 
nozzle and nuclear reactor core. This contribution can be 
diminished by creating an optically thick protective fluid 
layer enveloping the hot core gases as they expand through 
the nozzle. This protective fluid layer must exist adja- 
cent to the nozzle wall and may be formed by injection of 
a suitable cooling fluid through the wall. In addition 
to providing a radiation shield, the injected fluid must 
also diminish the convective heat transfer through the 
mechanism of transpiration cooling. 

To facilitate injection of the coolant fluid, the 
nozzle can be fabricated from a porous material. It is 
expected that injection on the order of 10 to 20 percent 
of the main stream mass flow may be required to provide 
corrective and radiation protection. In view of this 
large amount of injection, the changes in the flow-field 
characteristics due to injection need to be investigated. 

To the authors' knowledge, very little information is 
available on flow systems utilizing such high injection 
flow rates through a porous wall. The bulk of the exist- 
ing work applies to boundary layer flows in which the 
injection velocities are less than one percent of those in 
the external flow. These prior studies have been further 
limited to flows with small streamwise pressure gradients. 
Therefore, they do not apply directly to internal flows 
with large streamwise gradients in pressure and density, 
such as are found in supersonic nozzles. 

In the present report, an experimental approach is 
adopted for the study of nozzle flows with large wall- 


2 



transpiration rates. Based on the hydraulic analogy, a 
compressible gas flow is modeled with a water flow having 
a free surface. In its simplest form, the analogy applies 
between one-dimensional, open-channel liquid flow and 
isentropic, one-dimensional, internal gas flow. The flow 
variables in one-dimensional open-channel flow are local 
height and Froude number. They vary in the flow direction 
which is perpendicular to local gravitational acceleration. 
The local height and Froude number may be related to an 
equivalent local pressure and Mach number in an isentropic, 
one-dimensional nozzle gas flow. The analogy is exact for 
a perfect gas with a specific heat ratio y of 2. If the 
channel width-ratio (local width to throat width) is set 
equal to the nozzle area ratio (local area to throat area), 
the Froude number is equal to the Mach number, and the 
square of the height-ratio (local to stagnation) is equal 
to pressure ratio (local to stagnation). Once the Mach 
number or pressure ratio is found in this manner, the other 
properties may be found directly from the isentropic 
tables. For gases with different y, appropriate correc- 
tion factors exist (2). Previous reports in this series 
have described the analogy in more detail (3,4). 

The hydraulic analogy does not directly apply to a 
flow with injection. Nevertheless, it is useful in a qual- 
itative analysis of this type of flow. Such an analysis 
is performed here. It is simplified due to the assumption 
of inviscid and uniform flow. That is, the effects of 
viscosity and the flow variation across the nozzle are 
neglected. At any cross-section the flow has plug-flow 
characteristics with constant liquid height in the trans- 
verse direction. 

In the sections to follow, the emphasis is placed on 
studying the (a) Froude (Mach) number distribution, 

(b) shifting of the hydrodynamic throat (position where 
F and M = 1), and (c) thickness of the Injected fluid 
layer, all as functions of the injection rate. The calcu- 
lated results for one-dimensional flow are compared with 
experimental data obtained on the water table. Results 
presented here are for injection through a uniformly por- 
ous wall. Results for injection through a wall with dis- 
crete slots are presented in a subsequent report. 

Some results are also presented for injection of dif- 
ferent gases (air, CO 2 , Freon-12) into the transonic re- 
gion of a plane porous nozzle operating with air as the 
main flow. The qualitative agreement between the 
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incompressible water table results and those obtained in 
the wind tunnel is good. 

It should be re-emphasized that the primary objective 
of this work is to determine general changes in the flow- 
field characteristics caused by high rates of wall injec- 
tion, and not the precise prediction of the flow field. 
Therefore, the analytical results are used only as an aid 
in interpreting the experimental findings. 


HYDRAULIC ANALOGY 


The hydraulic analogy was first discussed in detail 
in Preiswerk (5,6). Loh (7) has reprinted much of this 
earlier work and has extended the analogy to the case of 
unsteady flow. In the analogy, an assumption is made that 
accelerations in the vertical direction are negligible. To 
approximately meet this assumption, it is necessary to 
operate the water table with modest stagnation heights in 
the upstream region where the velocity approaches zero. 

For the porous wall nozzle, this is difficult since one 
must have sufficient height in the supercri tical region 
after the throat for controlled blowing to be visualized. 
Thus wi th . f ejection , it is necessary to operate with stag- 
nation heights on the order of 2.5 to 3 inches instead of 
the more desirable 1 to 1.5 inches. 

To investigate the effect of stagnation height on the 
distribution of Froude number in a nozzle, a series of 
test were made with gradually decreasing stagnation heights 
These results are shown in Figure 1. 

The Froude number is analogous to the Mach number in 
compressible isentropic gas flows. These data were ob- 
tained without injection through the nozzle wall, the con- 
tour of which is shown in the figure. This nozzle has a 
45° inlet half angle and a 15° half angle in the down- 
stream section. 


As expected, the results with different stagnation 
heights become essentially indistinguishable as H 0 
approaches one inch. These results indicate that vertical 
acceleration effects are important for stagnation water 
deptns greater than approximately 1.5 inches. Preiswerk 
hai found similar results for an impermeable-walled nozzle. 


M?rh rvumbe. 


Figure 2 shows these same results plotted in terms of 
sir on- equivalent -gas wi Lh y - 1.4. Reference 
(4) has indicated the procedure by which the Froude number 
resul ts (Froude corresponds to the Mach number in a gas 
with y - 2 ) can be corrected to give Mach number results 
for other values of y. 


Briefly the correction follows from equating the area 
ratio formula for the compressible-gas flow to the width- 
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ratio result for water-table flow. These are: 


A - Ir / / 1 + X-l m 2 \ “i 

A th M u y+1 M 2 n u 


(y + 1 )/2(y“1 ) 


0 ) 


and 


fi. s lr( 1 /] + £ 2 \ *] 
B th F u 3 ' + 2 ' J 


3 / 2 


( 2 ) 


One can quickly see that in the analogy, P corresponds to 
M when y = 2 for the analogous gas. 

If one equates these two expressions, a single equa- 
tion results which can be solved using specified values of 
Y and F. The value of M so obtained is called the "cor- 
rected Froude number" and it is this value which is shown 
in Figure 2. Details for obtaining the correction are 
given in (4). Cuffel et. al . (8) have obtained Mach num- 
ber distributions in an axi symmetric nozzle having conver- 
gence and divergence half-angles of 45° and 15°, respec- 
tively. The centerline Froude numbers from the water 
table, corrected to y ~ 1*4, show excellent agreement 
with the Mach number results from (8) {see Reference 4, 
Figure ?). it is worth noting that the good agreement was 
obtained even with a stagnation height of 2.8 inches on 
the water table. Thus, we may conclude that for th case 
of no-injection, excellent agreement occurs between the 
incompressible open-channel flow and compressible flow in 
an axisymmetric nozzle. 

The question now arises as to whether or net the 
analogy can be applied to the nozzle flow with wall-injec- 
tion. To answer this, one must first obtain the governing 
equations for flow on the water-table and flow of a com- 
pressible gas. These equations are well known for the 
case of flow through impermeable-walled nozzles. However, 
they have yet fo be ved for the case where fluid is 
injected from the nozzle wall into the main stream flow. 
The derivation of the water-table equations will be pre- 
sented in this first section. The next section contains 
the derivation of the gas-flow equations, following which 
the basis for the analogy is discussed. 
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ANALYSIS OF ..WATER-TABLE FLOW WITH INJECTION 


Derivation of Equations 

The analytical approach taken here is similar to that 
adopted for the no-injection case. That is, the principles 
of mass continuity and flow momentum constitute the essen- 
tial basis for the analysis. As in the no-injection case, 
the character! Stic parameter for the open-channel liauid 
flow with injection is Froude number. However, an addi- 
tional parameter must be introduced for the injected fluid. 
This is the wall Froude number, F w , which is defined as 
the ratio of the injection velocity and (gH) 1/2 , where H 
is the local liquid height in the channel, and g is the 
local gravitational acceleration. Recall that Froude num- 
ber, F, is defined in the same manner, with the injection 
velocity replaced by local main fluid flow velocity. 

The injection is normal to the wall and encompasses 
the entire height of the area wetted by the main stream. 
That is, the injected fluid is always in contact with the 
main stream flow, and never above it. The flow is assumed 
to be steady, inviscid, incompressible, free of surface 
tension, with a hydrostatic pressure distribution and uni- 
form velocity at any cross section. For simplicity, the 
flow will be taken to be one-dimensional. 

Figure 3 depicts the flow boundaries and the elemen- 
tal control volume. The following terminology is used: 


V - flow velocity in x-direction 

v. - injection velocity 

p - fluid density 

P - hydrostatic pressure 

H - liquid height in the channel 

B - channel width 

0 - channel wall angle 
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The continuity equation can be expressed in the 
simple form, 


d(pVHB) - 


( 3 ) 


where the term on the left hand side is the mass flow- 
rate differential between downstream and upstream control 
surfaces in the main-stream fluid, while the right-hand 
side expresses the total flow rate injected into the con- 
trol volume (from both sides of the channel). In Equation 
(3) all of the variables, including v w , may be only a 
function of x, the flow direction. 

The contributing factors to the x-momentum are: 


i) Momentum of the main stream flui^ flowing 
into the control volume, namely mV or pV 2 HB. 

11) x-eomponent of the momentum of the injected 
fluid 

Hmre^w sin 8 - p w Vw 2 H dB 

where the bracketed term is the injection 
area per wall. 

iii) Momentum outflow, given by the first two 

terms of the Taylor series expansion about 
the x-location corresponding to the upstream 
control surface. This is 

P V 2 HB + d ( pV 2 H8 ) . 


Thus, the total increase in x-momentum of the flow coming 
out of the control volume, that is, the difference between 
the outflow and inflow momentunis, becomes 


d(pV z HB) - p w v w 2 H dB 


Since we have assumed inviscid flow, the net force acting 
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an the control volume is the change of the hydrostatic 
pressure force due to the fluid height change; thac is, 


- pgBH dH. 


This is a positive quantity since H decreases in the x- 
direction. Thus, the momentum equation in its entirety, 
becomes 


- pgBH dH = d(pV 2 HB) - p w v w 2 H dB 
This also may be written as 

- pgBH $ = pVHB £ ♦ V ^(pVHB) - P W * W »H $ 


From the continuity relation, the second term on the right 
hand side of the above equation is 


2p w v w VH/ c °s 8 


Also from continuity, 


dV 

dx 


n 


1 dH A l dB 
B dx 


r i un 

[ H HIT 


] + 


■...2p.w.yw... 

pBcos 8 


Here we have noted that dp/dx s 0 since the main stream 
fluid is taken to be incompressible. 

From the channel geometry, 


dB/dx » 2 tan 8 
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With the above identities, the momentum equation reduces 
to 


(1 


Us M 

gH' cx 



y-i_ + ^(.pM/-pivw5 

gH 


tan 0 


4 ( P W P ) Vw.V. - 1 
gH cos 0 J 


(A) 


As previously stated, the ratio of the liquid veloc- 
ity V to the propagation speed of a weak surface wave, 

/gH , is known as the Froude number, F. Flow is classified 
as subcritical, critical, or supercritical when F is re- 
spectively less than, equal to, or greater than unity. 

The ratio vJi/gTf is the wall Froude number and is given 
the symbol F w . 

Introducing F and F w into Equation (4) and noting 
that on the water table p w " p> one obtains 


dH _ 2_ 

dx ~ 1-F 2 


H 

Bcose 


C(F 2 


* *y) sin 9 - 2FF W ] 


(5) 


Equation (5) represents the variation of the liquid height 
in the channel including the effects of blowing. 

The variation of the Froude number with x location is 
found by combining the continuity relation written as 


HB 


dV 

dx 


+ VH 


dB + VB dH a 2y*H 
dx dx cos0 


with the expression 


dV _ V dF V_ dH. 

dx “ F dx 2H dx 


This is obtained by differentiating 
V 3 F/gff 
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The relation for dB/dx previously stated is also used. 

The final equation representing the Froude number distri- 
bution is 


dF 

dx 


3 F dH 2 

2 H lx B cos 0 


F sin 0 ) 


( 6 ) 


Equations (5) and (6) are coupled nonlinear first order 
differential equations. For a given nozzle geometry they 
can be integrated numerically to find the Froude number 
and liquid height variations as a function of x, the flow 
di rection . 

It is now of interest to examine the general flow 
characteristics revealed by the equations derived above. 
This will be taken up next, following which the corre- 
ponditg gas-flow equations will be derived. 


Flow Characteristics 

First consider the flow without injection, F w = 0. 
Equations (5) and (6) reduce to: 


dH 

dx 


2F 2 

T^T 2 ' 


g- tan e 


dx 


/ F 2 + 2 

(pcy 


) jj- tan 0 


(7) 

( 8 ) 


At critical flow (P = 1), the gradient in water height, 
dH/dx, must be continuous. This is analogous to the con- 
dition that the pressure gradient be continuous at M = 1 
in gas flows. Consequently, F = 1 must correspond to 
0=0, since otherwise the numerator of Equation (7) would 
not be zero. Recalling that 6 is the local slope of the 
nozzle contour (measured from the axis), it is seen that 
0=0 represents the position of the geometrical throat. 

It can be seen from Equation (8) that dF/dx is also con- 
tinuous at F = ( . 

Assume now that the subcritical region is upstream of 
the throat - , here by definition, F < 1 and 9 < 0 as indi- 
cated in Figure 3. From the Equations (7) and (8), in this 


11 



region dH/dx < 0 and dF/dx > 0. Consequently the flow 
must accelerate with diminishing height level toward the 
throat. Downstream of the throat, e > 0. If, in addi- 
tion, dH/dx < 0, then dF/dx > 0. On the other hand, if 
dH/dx > 0, then dF/dx < 0. Thus, for continually decreas- 
ing liquid height in the convergent-divergent channel, the 
flow must be accelerated to supercritical conditions 
(F > 1). We shall deal here with this category of flow. 

It is possible, however, to decelerate the flow after it 
has reached critical conditions at the throat. Thi". would 
be important in the analysis of diffusers. 

In a similar manner, we may arrive at the hydrodynamic 
throat conditions for the flew with injection. Here, the 
(shock free) critical flow position may be fouid from Equa- 
tion (5) assuming dH/dx to be finite for F = 1. Then the 
term in brackets must vanish. That is, when F = 1, 


d ♦ F„ tt «) sin 0 t - 2 f - 0 (9) 


Here F Wj t refers to the wall Froude number at the hydro- 
dynamic * throat , and 0t to the correspond! ng wall angle. 

The angle 8 t i s bounded by ±nr/2, and sin e t is required to 
be positive for Equation (9) to hold (i.e. F w# t cannot be 
negative). Thus, 0t must be positive. This means that 
the hydrodynami c throat (F = 1), for the case of wall 
injection, appears in the divergent portion of the nozzle, 
downstream of the geometric throat. 

The solution to Equation (9) is F w# ^ = (1 ± cos 0^)/ 
(sin 0^). For the nozzle under consideration the maximum 
value of is 15°. The corresponding values of F Wj t sat- 
isfying the Equation (9) are 0.131 and 7.59. The wall 
Froude number of 7.59 appears to be unrealistic since an 
injected fluid velocity of more than seven-fold the main 
stream velocity would totally disrupt the flow. Further- 
more, as 0-t approaches zero, the upper values of the wall 
Froude number get arbitrarily large. Thus, only the lower 
values of the wall Froude number, given by 


P _ 1 -cos Bf 
r w,t“ sin e t 


( 10 ) 


appear to have physical significance. 
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At F w equal to F w ^ the expression for ( d H/ dx ) t be- 
comes i ndetenni nant . however, its value may be determined 
with the aid of l'H6pital‘s rule. One finds 


dH) 2Hj_(2 F ( d_Fy d x ) + 2F w Ld,Fw/dx ) ) si n ej, + 

W t - 2FB cos 0 (dF/dx) 

8-*0 £ 

F-*l 


2 H { ( F 2 + Fw 2 )cos 6(d6/dx)} 
-2FB cos 0 (dF/dx) 

2H { 2 [ F ( d Fyj/ dx ) + F w ( d F / d x ) j } 
- 2FB cos 6 (aF/'dx) 


If the wall injection velocity is constant with x, then 


d£w = . 1 f.w dH 
cTx 7 H dx 


When this is inserted into Equation (11) together with 
Equation (6) defining dF/dx, one obtains a quadratic equa- 
tion in dH/dx. The final solution, after taking the limit 
and considerable algebra, is 


37) r ■ 3B33 TT { < 6F w - 6 si " e ' V s1 " e > 1 

[ ( 6 F. ( - 5 sin 9 - F,, 2 sin e) 2 + 
w w 

6 ( { 1 + F 2 ) B cos 2 0 (de/dx) - 
w 

4(F W - sin 0) 2 )] l/2 } 

with all variables being evaluated at the hydrodynamic 
throat . 


In the throat region the channel wall has a circular 



shape of radius R. Thus 


49. a . -1. _ 

dx R cos 0 


For F w,t given by Equation (10) the expression for the 
slope of the 1 i qu id-he i ght further reduces to 



4H 




- cos 


B sin 


{1 ± 


ri + 3/ B 

L 3*'R cos 0 (1 -cos 


- 2 )] 1 / 2 } 


( 12 ) 


The positive sign will be chosen io render the value 
(dH/dx)^ < 0, since for small 0 and B of the same order as 
R the second term under the radical sign is positive. 

Equations (5), (6), and (12) together with the known 
nozzle geometry yield sufficient information for numerical 
studies of the Froude number and the liquid height varia- 
tions along the nozzle. In addition, the thickness of the 
injected layer measured normal to the wall, 6 , may be 
estimated using the simple expression 


fi (x) 


total mass injected upstream of the x 

posi tion 

“ " ( 2 pH V/ cos 0) 


(13) 


This assumes that the injected fluid displaces the main 
flow adjacent tc the wall without mixing with the free- 
stream flow. This expression follows from the continuity 
equation for the incompressible injected fluid. In keep- 
ing with the one-dimensional idealization, the fluid in 
the injected layer is assumed to have the same velocity, V, 
as that in the main flow. 

Results of calculations using these equations will be 
presented shortly. 
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ANALYSIS OF COMPRESSIBLE-GAS FLOW WITH INJECTION 


Derivation of Equations 

Let us now look at the problem of a plane porous noz- 
zle with a compressible fluid injected through the wall 
into a compressible main scream fluid. 

Our model will assume the following: 


1 . inviscid flow. 

2. injection normal to the nozzle wall. 

3. perfect gas behavior (or Pv - 2RT can be 
used ) . 

4. a plane porous nozzle with injection only 
through two opposite sides. 

5. unit depth normal to the plane of flow, i.e. 
H = 1 . 


The continuity equation then becomes (see Figure 3 
for the geometry) 


|^ = ^(PVB) (14) 


The momentum equation can be developed in a manner similar 
to the incompressible case. We have 


£ F = Rate of Creation of Momentum 

A 


In the absence of shear: 


l F, 


1 5 


B dP 



‘he momentum flow into the control volume is pV 2 B 
plus the X-momentum brought in with the injected fluid. 
This latter quantity is 


° <«Vw dx/cos 0 > v „ si " e 


Since (dx/cos 6 ) (sin 6) - dB/2, we have 


(in inj)in = Vw 2 dB/2 


The complete momentum and force balance then becomes 
(injection through two sides) 

- B dP = ^ (pV 2 B)dx - p w v w 2 dS 


In differential form we have (note dB - 2 dx tan 0 ) 


dJL . tan e - 1 ( P VB) ^ ^ CoVB) (15) 


Substituting for the last term from the continuity equa- 
tion we find: 


ILL a ( v + an Q _ — V — ) _ a \j dV ( 16 ) 

dx B ' v w r n 0 cos 0 ^ pv 3x ‘ 


The energy equation can be developed in a similar manner. 
Since 


E . = E 

in out 
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and 


E,„ ■ pV8(V a /2 + h) + Z g^ w j| X (» w 2 /2 + h w ) 
E out “ pVB(V 2 /2 + h) + ^[pVB(V 2 /2 + h)J dx 


one obtains, 

|gff[(v w a /2) + h„] - ^[pVB(V a /2 ♦ h)3 (17) 


We can now manipulate these equations, using also the 
perfect gas equation, to find a relation for dT/dx. In a 
manner similar to that used for Equation (5) the effect of 
injection on the location of the hydrodynamic throat can 
then be determined. 

Combining the energy and continuity equations and 
assuming that Cp iW a c p we have 

♦ c (T w -T» - pVBCV $ ♦ c jjl] 

(18) 


From continuity 


0 B 37 + P v 37 + VB & = < ,9 > 


The perfect gas relation can be differentiated to give 



( 20 ) 


Substituting the momentum equation for dP/dx, one finds 
dp/dx in terms of dT/dx and dV/dx. This expression is 
then substituted into Equation (19) to find dV/dx. Noting 
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that dB/dx ■ 2 tan 9 we have 


dV . _jp»v M 
dx pBcos 0 


2Vtan 9 

§ 


v p dT 
p T dx " 


Jt 9 - */ cos 9 3 * RT 37 


or 


P-[] - Jjl] , -p*** . 2Vt an 9 V dT 

37 u RT J pbcos e 6 T 3x " 

v^tan 9 - V/cos 9] (21 ) 

Now substituting this expression into the energy equation. 
Equation (18), we find 

h7 [cVBc p * T (f-vvat) J ■ * «p ( V T)] 

' T^TTm ^ 7*5 - 


v w tan 0 


V/cos 0 ] 

( 22 ) 


Using 


M 2 - V 2 / y RT 
and 

V “ V/yRT 

(note is defined on T, not T^) we can reduce this 
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equation to some extent. For example, the bracketed term 
on the left hand side becomes 



Using this, an intermediate result is 


& ■ g rffiV - M2 > + 

Tc p(f w - , )J-ptee TT-V7'K1T (1 + yM “ ' 

yMM w sin 0) + ■ B^ f - ~ V ' VR ~ TT } (23) 


Further manipulation results in the final equation (note 
YR/c p ■ (r-1)): 


dT 

37 = 


II.-C-j.j J, U-Yl 2 lflwVw r MV-M* 
U-M 1 ) B 1 pVcos 6 L 2 



+ yM 2 - yMM w sin 


+ M 2 tan 6} 

(24) 


This rather complicated expression is the counterpart to 
Equation (5) obtained for the water-table flow. An expres 
sion for dM/dx could also be obtained with steps similar 
to those used to obtain Equation (6). 


Flow Characteristics 

At the hydrodynamic throat, M = 1. In order for 
dT/dx to be continuous there, the term in braces must 
approach a limit of zero. We then have, with M = 1, 


si " e t - c w e t> + 

t . i . T y/T-l) » M wt (l-y)] - 


1 9 


0 



Expanding, one obtains the cubic equation, 


M 3 t + 'T 1 f sin “ tfpH + 
w , t 1 -y W , t t ‘•(I -y) (1~Y ) w , 


2o 

+ iTTi^rr 


sin 0 , = 0 
t 


(25) 


Solving the cubic for as a function of 0 t wi11 reveal 

the location where M - 1 in the nozzle. A similar situa- 
tion arose in connection with F w ,t and 0t when F ~ 1. 
However, in that case a quadratic equation resulted for 

F w,t* 


Solution to the Cubic Equation 
Consider the equation 


x 3 H| x 2 t a 2 x + a 3 0 


If we let 


Q 3 lsLLjj_§«L Z 


9a ia»-27a^-2ai ? 
54 


then if D = Q 3 + R 2 , we have the following possibilities: 


(a) one real root and two complex conjugate 
roots if D > 0. 

(b) all real roots (and at least two are equal) 
if D = 0. 

(c) all roots real and unequal if D < 0. 


From our equation 


ai s sin e t 
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a 2 

a 3 


. | ; (v + 3)' f 2(Tw/T-l ) j 
Zsinilt p 

1“Y P W 


It becpmes algebraically cumbersome to develop expres 
sions for Q 3 and R , so we shall look at ai, a 2 » and a 3 , 
directly. Since for compressible fluid, y > 1, and wo 
are Interested in 0 £ 6 t < v/2, we have at < 0. For as, 


Note that 
2 ( T w / T - 1 


all cases of interest will have Tw/T <_ 

( y+3 ) > 0, and in particular |(y+3)j > 

Since (1-y) < 0, the numerator is positive and the denom- 
inator is negative. Therefore, a 2 > 0. Concerning a 3 , 
we see that a 3 < 0. 


In summary, at and a 3 are negative and a 2 is posi- 
tive. Since Q = (3a 2 -a£ 2 )/9 it can be either positive or 
negative. Therefore we must look at the magnitudes of the 
terms to establish the sign of Q and also of R. 

Let 8t" 10°, then for y * 1 .4 we have 


a, = sin 6 t = -(7) (0.173) = -1.2 


For T w << T we have 

a* ■ ‘ " 6 


So in order of magnitude 

g e 3 . a .2. -a . i . 2 . « 1 (is - i .4) a 1 .85 


or Q > 0 . 

Now, the sign of R can be established by looking at 
a 3 . We have 


a 3 


2fLn«t£_ < 

(l-y) Pw 


o 
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For a perfect gas p/p w c T^/Tf So for an actual rocket 
nozzle p/p w < 1. We shall pick, rather arbitrarily, 
p/pw u 1/4. Then 


a 3 ° ^-|47s 7^ (-25) « - 0.173 


and 

R o 9(-l . 2 ) ( 6 ) -27 ( - . 1 7 ) -2 ( -1 .2? 
54 

or 


R = 


-56.4 
54 “ “ 


Thus R < 0. Careful consideration of the possible range 
that the variables can take indicates that it is gener- 
ally true that R < 0. Recalling that 0 = Q 3 + R 2 , and 
since Q 3 > 0 and R 2 > 0, we have D > 0. Therefore, one 
root is real and two are complex conjugates. OetermTnTng 
the one real root will establish the location for M - 1. 

Because there are a large number of parameters in- 
volved, the solution for Mw 9 t will be obtained for select- 
ed values for T/T w , p/pw and y. The equation will be 
solved for the one rea root to determine the K w , t value 
as a function of B^. 

The solution is presented in Table I. The geometry 
is noted in Figure 4. The procedure Involves picking 
vaiues of sin 0^ and obtaining M Wj t by successive trials. 
Thus for T w /T s p/pw = 1*0 and y - 1.4, we find that a 
value of M w = 0.090 will produce M = 1 at 0^ s 11 ° 2 0 ' . 

This angle is measured counterclockwise from the geo- 
metric throat (where 8 = 0). The effect of injection is 
therefore to cause the hydrodynamic throat to move down- 
stream from the location of A* associated with the geo- 
metry, (i.e. for non-porous nozzles, M = 1 when A = A*). 


*This assumes that the injectant gas is the same as the 
propellant gas, for example hydrogen, and enters the 
nozzle at the local static pressure. 
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TABLE I LOCATION OF M “ 1 WITH INJECTION 





q 

t 


sin o 

T /T 
w 

“ 1 - p/p w u 

1.0, y 

s 1.4. 




0.045 


5“ 45' 


0.100 


0.090 


11° ; 

20' 


0.200 


0.130 


17° 

28' 


0.300 

V T 

s 1 / 2 , p/p w a 

1/2, y 

° 1.' 

4. 




0.0294 


5° 45' 


0.100 


0.0588 


ir 

20' 


0.200 


0.0882 


17° 

28' 


0.300 

V T 

= 0.1, p/p w » 

0.1, Y 

= 1. 

4. 




0.0077 


5° 

45' 


0.100 


0.01 54 


11° 

20' 


0.200 


0.0231 


1 7° 

28' 


0.300 

V T 

53 0.01 , p/p w : 

3 0.01 , 

Y “ 

1 .4 

* 



0.00083 


5° 

45' 


0.100 


0.00166 


11° 

20' 


0.200 


0 . 00249 


17° 

28' 


0.300 

■HI 

£ 

— ! 

s o, p/p w s 0 

, Y s 1 

. 4 , resul t is: 



M., «* 0 for 
w 

all e t 

• 




(Effect of y.) 






T /T 
w' 

s 1.0, p/p w = 

1.0. 






y M 

r ' w 


0 t 


M w 

0 t 


1.2 0.047 

5 e 

45' 


0.142 

17° 28' 


1.4 0.045 

5 e 

45' 


0.130 

17° 28' 


2.0 0.040 

5° 

45* 


0.120 

17° 28' 
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This flow characteristic is in qualitative agreement 
with tnat revealed by the analysis of the porous-nozzle 
flow on the water table. Thus, one can conclude that 
incompressible flow studies on the water table can be 
used to investigate (qualitatively) the effects of mas- 
sive injection on the dynamics of compressible flow in 
nozzles. 

Before discussing the flow analogy further, we can 
complete this development by applying these results to a 
hydrogen-propell ed gas-core nuclear rocket. The inlet 
pressure to the nozzle will be taken as 1000 atm (Kascak, 
(9); Ragsdale and Wiilis (1)); the gas inlet temperature 
is 20 , 000° R (11 , 1 00°K) . The mass flow rate is 10 lbm/sec 
(Taylor et. al . (10)). Properties of hydrogen are taken 
from Kubin and Presley (11). 

At the throat P/P 0 - 0.5 and !/T 0 * 0.8 assuming 
isentropic flow. So we estimate that 


T* - 1 6 , 000° R 
P* s 500 atm 

2 a 1 .90 (pv = 2RT ) 

p - 0.0448 1 bm/f t 3 

Injected coolant is assumed to be hydrogen at 500 atm 
and 540°R. This pressure is picked so that the wall pres- 
sure approximates the main stream pressure. The temper- 
ature represents a nominal value. Under these conditions 
p s 2.56 lbm/ft 3 . Therefore 

T w /T = 0.034 
p/p w = 0,0175 

Values of y are 1.40 at the wall and 1.50 in the main 
stream. Taking y - 1 .5 and using the above data, we can 
solve the cubic equation for M w> t to find the location 
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fo r M = 1 . We f i nd 


5 e 45' 0.00136 

17 s 28' 0.0040 


Thus, a wall (injection) Mach number of 0.00136 will move 
the hydrodynamic throat 5° 45' downstream from the A* 
position. The blowing parameter, X, is given as 
X = (pwV w )/pV. At the throat, for 0.001 36, we have 


X = (1/0.01 75) (0.001 36) = 0.078. 


The above calculations illustrate the manner in which the 
throat characteristics can be estimated for an actual 
nozzle with wall injection. 
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DISCUSSION QF HYDRAULIC ANALOGY IN POROUS NOZZLES 


The basis for the hydraulic analogy, if one exists, 
is obtained through a comparison of the flow equations 
for the water-table and compressible gas flows. For our 
purposes, Equations (5) and (24) are sufficient. 

It is of interest to first specialize these equa- 
tions for the case of an impermeable nozzle (M w = F w = 0). 
One obtains for the water-table and gas flows, respec- 
tively , 


dH _ 2F 2 H . 

37 - TT* B tan 6 


(26) 


dT_ 

dx 



0 


(27) 


Obviously, the substitutions M - F, y = 2, and T = H in 
Equation (27) produce Equation (26) exactly. On this 
basis we may say that an exact analogy exists when y ~ 2 
and the nozzle wall is non-porous. This fact has been 
known for decades and is called the hydraulic analogy. 

For the case of injection, Equations (5) and (24) 
must be compared with finite M w and F w . While there are 
similarities between the two equations, there are no 
simple substitutions which reduce Equation ( 24 ) to Equa- 
tion (5). For this reason, a direct analogy does not 
appear to exist between compressible nozzle flow and incom- 
pressible open channel flows with wall injection. This 
may be due to the fact that with injection, the flow does 
not follow an isentropic process except in the trivial 
case of injection of the same fluid as the main stream with 
v w = V , T w = T, p w = p, P w = P. Also, the injection must 
be parallel to the main stream. 

For injection through a porous wall the aforementioned 
conditions cannot be met. For injection through a slot, 
however, the conditions can be approached. Nevertheless, 
even then one would presumably not inject fluid at the same 
temperature as the main stream. 
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In conclusion, the nozzle flow is not isentropic when 
there is uniform wall transpiration, and an exact analogy 
does not hold. However this does not mean that analogous 
behavior does not occur. In fact, such behavior does occur 
and the required equations to investigate this have already 
been developed. The equations show that for compressible 
flow, the effect of injection is to move the position of 
the sonic point downstream of the geometric throat. The 
same effects occur for the incompressible open channel 
flow. 


Remits from the experimental portion of this inves- 
tigation are presented and discussed in the next section. 
They show that wall transpiration indeed shifts the posi- 
tion at which F = 1 in the downstream direction. This is 
in agreement with the trend already found from the one- 
dimensional flow equations. Because of the close corre- 
spondence between the experimental findings and the one- 
dimensional flow analysis, complete numerical solutions 
have been obtained for the water-table equations, Equations 
(5) and (6). These are presented and discussed in a 
later section. 
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EXPERIMENTAL RESULTS WATER TABLE 


Qual i tati ve 

Because the emphasis in the experimental work was on 
observations of flow phenomena, a large number of colored 
slides, black and white stills, and a lesser quantity of 
16 mm movies were obtained for flow with injection. Re- 
production cost limit the presentation of these results to 
black and white stills. Figures 6, 7, 8, 9, 10 present 
results for the effect of injection on the flow field in 
the nozzle. The conditions under which the figures were 
obtained are listed in Table II. Figure 5 shows the noz- 
zle as positioned on the water table. Blowing rates can 
be individually controlled in the three sections labeled 
I, II, and III on the Figure. 

The flow visualization technique, described in detail 
in (4), produces color differences between fluid streams 
without the use of permanent dyes. This is accomplished 
through pH control of the two streams and a suitable acid- 
base indicator ( bromothymol -bl ue ) mixed in the fluid. The 
injected fluid is made basic with the addition of sodium 
hydroxide; the main flow is an aqueous solution of acetic 
acid. 


Figure 6 notes flow without injection. The expansion 
fan is evident in the divergent portion of the nozzle. 

The variation in water height (in the direction perpen- 
dicular to the plane) produces a pronounced distortion of 
tf‘ i straight lines ruled at one inch intervals along the 
t e surface. (The lines ruled across the channel are 
not at one inch intervals). 

Figures 7a and 7b compare injection at low and high 
rates in section I. Since here the main stream velocity 
is low the injected fluid can penetrate well into the main 
flow. Thus injection of 1.85% (per wall) of the inlet 
mass flow (Figure 7a) produces a thick, well defined layer 
in the converging portion of the nozzle. 

Figures 8a and 8b compare injection at low and high 
rates in the second section of the nozzle. With 1.45% 
injection per wall only a very thin protective layer is 
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TABLE II FLOW CONDITIONS FOR FIGURES 6 TO 10 



Injection Conditions 


4 * 

Figure 

Section I 
"i^lnlet* 

Section II 
'"il'Snlet 

Section III 
m III /m inlet 

6 

0 

0 

0 

7a 

3.9% 

0 

0 

7b 

8 . 2 % 

0 

0 

8a 

0 

2.9% 

0 

8b 

0 

7.6% 

0 

9a 

3.1% 

5.7% 

0 

9b 

3.1% 

3.1% 

3.1% 

10a 

1 .5% 

1 .5% 

1 .5% 

10b 

6.6% 

6.6% 

3.3% 


+ Plus-x Film, ASA 125, top and bottom lighting, 1/125 at 
f8 exposure. 

*Injection from both sides so 8.2% represents 4.1% from 
each side. 
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provided at the throat. With substantial injection of 
3.81 of the inlet mass flow, Figure 8b, a well defined 
thick layer is provided. In this region the main stream 
velocity is about 1.6 ft/sec at the throat, and the 
injected flow does not penetrate very far into the main 
fl ow. 


Figures 9a and 9b contrast blowing in the first two 
sections with blowing in all three sections. Blowing in 
the third section produces pronounced enhancement of the 
surface waves in the divergent portion of the channel; 
compare Figure 6 with Figure 9b. 

Figures 9b, 10a, and 10b are for blowing in all three 
sections at total injection mass flows of 9 . 3 %, 4.5%, and 
16.5%, respectively. The injected layer maintains its 
integrity rather well considering the effects of vertical 
accelerations disturbing the flow in the transonic region 
of the nozzle. 

Even massive injection does not adversely affect the 
main stream flow. The injected layer remains intact, and 
while mixing at the interface occurs, the main stream does 
not penetrate to the wall in any of the cases with injec- 
tion in all three sections. It will be seen that these 
pictures correspond closely to those obtained with the 
gas flow test section. 

It seems appropriate at this point to comment on the 
magnitude of the mass injection. The usual parameter to 
describe injection is 


X = 


o v / pV 
y w w 


This is suitable for external flow but for internal flow 
a better parameter might be 


X 


int 




to 


If p w and v w are functions of position, then an integrated 
value is required along the wall. The reason that X is 
not fully descriptive for internal flow is that a small 


30 



value can still lead to large ratios of if the 

injection occurs over a big area. Thus a more appropriate 
quantity must include the area for wall injection. 

Since water table results are for incompressible flow, 
p w s p and X reduces to 


X 


v K /v = rjr 


Typically at the throat Vt =1.6 ft/sec and a typical 
value for the wall velocity is v w s 0.03 ft/sec. Thus 
X = 0.03/1.6 * 0.02. 


Quanti tati ve 

Figure 11 shows the no-injection Froude number as a 
function of position in the nozzle. The two-dimension- 
ality of the flow field is evident. Froude number distri- 
butions are obtained by measuring the local water depth 
to obtain /gTJ". The water depth is determined by measuring 
with a pointed depth micrometer the distance to the water 
surface and the distance to the table surface. The dif- 
ference is the local water height. The Froude number is 
then determined from Equation (28). 

Also shown in Figure 11 is the corrected Froude 
number (corrected to a Mach number corresponding to 
Y - 1.4). The correction, previously discussed, is partic- 
ularly significant at high values of F. 

Figure 12 shows Froude distributions with and without 
blowing. The solid curves correspond to the case of no 
injection, and the dashed curves show the effect of wall 
blowing. The blowing conditions are similar to those in 
Figure 10b and the relative rates of injection through the 
three sections are as indicated. It is seen that the 
effect of injection is to move the critical "line" F - 1 
downstream. 
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INCOMPRESSIBLE SOLUTION 


COMPUTATIONAL PROCEDURE AND RESULTS 


The governing equations for incompressible flow with 
injection were derived previously. Equations (5) and (6} 
are coupled nonlinear, first order, differential equations 
which can be solved to yield numerical predictions for 
Froude number and the liquid-height spatial distributions. 

Runge-Kutta 1 s fourth-order integration scheme was used 
to integrate these nonlinear equations. The error intro- 
duced is known to be of order (Ax) s . The Ax increment 
was varied from 0.0001 to 0.01 with the minimum value 
being at the throat region where the slope changes are the 
largest. The curves F(x) and H(x) are very smooth and 
hence the error introduced is not significant. 

To insure the continuity of dH/dx and dF/dx at the 
hydrodynarni c throat, integration commenced at that point. 

By picking the value of 0^ between 0° and 15°, the value 
of F w t was computed from Equation (1 0) and then the 
injection velocity from v w = F W} t^gHt* The value of 
liquid height at the throat was chosen to correspond 
closely to the experimental heights of order 1 to 1.5 
inches. However, it has no bearing on the dimensionless 
height ratio (local height to stagnation height) or the 
local Froude number. For the given 0t, the height <.nd 
Froude number slopes were determined from Equations (12) 
and (6), respectively, and the values of H and F calculated 
a Ax increment away. This was done for several minimum 
step sizes, sufficiently away from the region of insta- 
bility of Equation (5). Afterwards, Equations (5) and (6) 
were exclusively used in the Runge-Kutta scheme of inte- 
gration, first upstream of the hydrodynamic throat and 
then downstream of it. The injection velocity (and thus 
F w ) could be varied at each step (or set equal to zero, if 
so desired, as was done for the region upstream of the 
straight 45° convergence portion of the nozzle). In this 
manner, comparisons could be made with experimental data. 

The calculations were not carried to the extreme 
upstream region. Rather the stagnation height was 
determined from the expression 
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(28) 


H 0 s H (x ) [1 + 0.5 F 2 (x) ] 


which is a statement of the conservation of energy along 
the streamline. The numerical results are presented 
graphically. The nozzle of interest, already noted in 
Figure 5, is divided into three sections in which blowing 
rates can be individually controlled. These are labeled 
I, II, III in the Figures. 

Mass injection has a significant effect on the Froude 
number, especially in the supercritical reglon.as may be 
seen.in Figure 13. The parameters are St and minj/roi nl et ■ 
For Hiinj/ m inlet> a value of 0.13 corresponds to total 
injectfon equal to 13% of the inlet mass flow. Since the 
injected mass x-momentum is onTy~a — smal 1 fraction of the 
main stream momentum (in the convergent nozzle section the 
injected x-component is negative), the mainstream mass 
must impart some of its kinetic energy to accelerate the 
injected mass. This manifests itself in a lower Froude 
number as compared to flow without injection. Since the 
kinetic energy is proportional to the square of the veloc- 
ity, the energy expended by the main stream mass to accel- 
erate the injected mass in the supercritical region is 
accordingly much greater than in the subcritical region. 
Hence divergence of the Froude lines from the no-injection 
case is more pronounced when Injection takes place in the 
supercritical region. Note that the effect of Injection 
is to move the position at which F - 1 downstream of the 
geometric throat. 

Figure 13 also depicts the H/H 0 curves for two injec- 
tion ratios, again defined as a percentage of the main- 
stream inlet flow rate. The increase in height ratio is 
confined mainly to the supercri tical region. The effect 
of injection is to decrease the mainstream velocity at a 
given x-location. This is accompanied by a local increase 
in the liquid height since continuity must be preserved. 
Thus the effect on the Froude number is compounded, because 
F is proportional to V and inversely proportional to /FT. 
This is clearly exhibited in the figure, the influence on 
H/H 0 being much smaller than on F. 

Figure 14 indicates results for F for the case of 
injection in sections I and II of the porous nozzle. For 
a total injection equal to 31£ of the inlet flow, the 
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position at which F - 1 is shifted downstream about 0.25 
Inches. This corresponds to the "sonic" line occurring at 
0 - 5®, as opposed to 8 s 0°. 

Figure 15 indicates the blown-layer thickness, A, 
as a function of position for various flow conditions. 

These results were calculated from Equation (13) and 
assume no mixing with the free stream flow. They there- 
fore represent minimum distances for which the effects of 
blowing are felt. 

The conditions for which these results were obtained 
correspond to blowing in sections I, II, and III , with equal 
v w in all three sections. In particular, for a total 
Injection of 1 3 % of the inlet mass flow, the injection 
velocity was 0.34 inches/sec. For an injection of 20%, 
v w a 0.51 in/sec and for 39%, v w - 0.60 in/sec. Since the 
water height is markedly dropping through the nozzle, the 
injected mass flow is not equal in each section. The shape 
of the o ( x ) curves is as noted because of this change in 
water height. In a gas nozzle, if injection occurred at 
the local pressure and with the same wall velocity at all 
locations in the nozzle, a similar shape should result. 

This is because the fluid density decreases in the flow 
direction. 

Figure 12 has already presented experimental results 
for the Froude number distribution with and without blow- 
ing. The data on this figure are for a total injection 
rate of 21.1% of the inlet flow, with an injection equal 
to 16.6% of the inlet flow in the first two sections. 

The F - 1 position along the centerline is shifted about 
0.2 inches downstream. The results of Figure 14 indicate 
a similar downstream shift in the "sonic" position. Thus 
the experiments and analysis show satisfactory agreement 
for the incompressi ble case with injection. 
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EXPERIMENTAL RESULTS - GAS TUNNEL 


Late in the experimental program it was decided to 
utilize the Aerospace and Mechanical Engineering Depart- 
ment's small supersonic facility for a flow visualization 
study of gas injection into a main stream air flow. 

Schlieren pictures were obtained with a spark light source 
for CO 2 and Freon~i2 injected through the porous walls of 
a nozzle geometrically similar to that used on the water 
table. Pictures for injection of air were obtained using 
a continuous light source. 

The wind tunnel is an in-draft type facility, and 
about 45 seconds of run time can be obtained. No quanti- 
tative results were attempted for this report. However a 
subsequent report on a nozzle with discrete slot injection 
includes some transient pressure results. In this respect 
the water table offers a real advantage since detailed fluid 
studies can be made over several hours of operation at the 
same flow conditions. 

Figure 16 shows the porous-walled nozzle used in the 
supersonic tunnel. Blowing rates could be independently 
controlled in each of the three sections. The section 
dividers are shown cross-hatched. 

Table III summarizes the flow conditions for Figures 
17, 18, 19, 20, and 21. Figure 17 shows the flow with no 
injection. Figure 18 is a composite picture corresponding 
to injection of air into air. With air as the injectant, 
only small density gradients occur. Thus the injected 
layer is visible only in the supersonic portion of the noz- 
zle. The upper half of the figure shows details of the 
flow with no wall injection. The lower half corresponds to 
wall blowing with a total injection of 16.5% of the inlet 
mass flow, distributed as 3.3% for each side of section I, 
3.3% for each side of section II and, 1.65% for each side 
of section III. The oblique wave pattern in the divergent 
portion of the nozzle is Qualitatively the same as that 
found on the water table (Figure 7). 

Figures 19a and b show Freon-12 injection through all 
three sections. Figure 19a shows the effect of a low in- 
jection rate and 19b a high rate. Figures 20a, b, and c 
are for COg injection. Figure 20a has an injection of 
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TABLE III FLOW CONDITIONS FOR FIGURES 17 TO 21 


Injection Conditions 


Figure 

Section I* 

Section II 

Section III 

l "l /,,1 inlet + 

m II /m inlet 

m 1 1 1 / ,n i n i e t 

1 7 

02 

02 

OS 

18 

6.6 

6.6 

3.3 

(air) 

19a 

3 

3 

3 

{ freon) 

19b 

12 

12 

1 2 

20a 

3 

3 

3 

(co 2 ) 

20b 

6 

6 

6 

20c 

12.5 

0 

0 

21a 

12 

0 

0 

(freon) 

21 b 

0 

12 

0 

21c 

12 

12 

0 

* injection 

from both 

sides so 6% is 3 % 

from each side. 

+ cal cul ated 

from one- 

dimensional flow 

theory . 
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1 , B% of the Inlet mass flow (calculated from one-dimen- 
sional flow theory without blowing) in each section and on 
each side. Thus the overall injection rate is 9% of the 
inlet mass flow. Figure 20b shows COg injection at 3% of 
the inlet mass flow in each section and both sides (18% 
total injection). Figure 20c shows an injection from sec- 
tion I only. The injection rate is 12.5% (total from both 
sides) of the inlet flow. Figures 21a, b, and c are for 
Freon-12 injection from sections I, II, and I and II, 
respectively. 

For all of these cases it can be seen that injection 
does not disturb the mainstream flow pattern markedly from 
that found with no injection. 

It is possible to calculate an approximate value for 
the blowing parameter, X, from the data. The mass flow to 
each section was metered. Knowing the injection area, the 
product ( o v ) j n j can thus be calculated. The product (pV)^ 
can be calculated from one dimensional theory without in- 
jection. These have been obtained using the area ratios 
at the midpoint of the injection sections. The values of 
X - ( pv) n j/ ( pVju, thus obtained are given in Table IV. 

For section I, (pV) w is quite low so A is very large. 

In turn, this means that a thick injected layer can be 

expected. This is verified in Figure 20c. This layer is 
carried smoothly through the throat, and with some addi- 
tional, but lesser, injection in the throat region, should 
provide good thermal protection for the wall. Not sur- 
prisingly, X is low in the throat region, and injection 
there provides only a very thin layer; see Figure 21b. 

The qualitative agreement between these figures for 
compressible flow and the previous figures for flow on the 

water table is good. On the water table, injection in the 

subsonic region also produced a thick layer, but transonic 
injection results in only a thin layer. 

It is important to note that the thickness of the 
injected layer is not a linear function of the mass injec- 
ted. The low injection rates produce almost as thick a 
layer as do the very large injection rates (i.e. Figure 
20a with 9% compared to Figure 20b with 18%). This is 
probably due to the increase in the main stream pressure 
(at a given position) as the amount of injection increases. 
Thus doubling the injected mass flow also increases the 


37 



TABU IV BLOWING PARAMETER, X, FOR SEVERAL 
MASS INJECTION RATIOS 

A = (pv) . n j/ ( PV)^ 


m 1 n j 
"'inlet 


Section 

1% 

3% 

6% 

10% 

I 

0.095 

0.28 

0.56 

0.95 

II 

0.0094 

0.028 

0.056 

0.094 

III 

0.0115 

0.0345 

0.069 

0.115 


38 



density (and thus the momentum for a given velocity) of 
the mainstream fluid. As a result, X is not doubled. Thi 
parameter is believed to be more important in correlating 
penetration depth than the mass injection ratio. 



DISCUSSION OF RESULTS 


The results of the one-diniensional analysis for both 
incompressible and compressible flow have indicated that 
the position of the hydrodynamic throat is moved down- 
stream of the geometric throat because of injection. The 
more the injection, the greater the change in the position 
of the hydrodynamic throat. Without injection it is well 
known that the two-dimensionality of compressible flow in 
a nozzle causes the sonic position to be downstream of the 
minimum flow area. The effect of injection is to exagger- 
ate this effect. 

While no direct analogy could be established between 
compressible and incompressible flow, both with injection, 
analogous behavior was noted. For the case of no injec- 
tion the hydraulic analogy indicates that Froude number 
results are equivalent to Mach number results for a gas 
with y * 2. With injection such a simple relationship 
could not be established, but the results, both analytical 
and experimental, clearly indicate that water table 
results can be used to predict qualitatively what changes 
will take place in a gas nozzle with injection. 

The incompressible experimental results for the.porous 
nozzle show that even with massive injection {where mtnj 1s 
of the order of 20% of the inlet mass flow) the main stream 
flow is not seriously disturbed. A one-dimensional analysis 
could probably be used to estimate the flow field if the 
area used in the calculation were the actual geometric area 
less the area required for the injected layer. 

The compressible experimental results, using a 
Schlieren system to visualize the flow field, also indi- 
cate that massive blowing does not cause unusual phenom- 
ena to occur in the flow. The blown layer stays close 
to the wall and does not, for example, mix rapidly with 
the core flow. A potential problem was, however, noted 
which affects the thermal protection at the nozzle throat. 

It is difficult with a porous nozzle to obtain a thick 
blown layer at the throat region. The extreme axial 
acceleration at this poin 4- causes injected fluid to be 
moved downstream before penetrating very far into the 
mainstream. To obtain a thick layer at the throat, it is 
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necessary to have large injection upstream of the throat 
so that the blown layer can be swept through the throat 
intact and adjacent to the nozzle wall. Alternatively, 
one might obtain a thick wall layer at the throat by 
injecting fluid through a single (or multiple} slot 
upstream of the throat. 


CONCLUSIONS 


1. The Froude number at a given location in the super- 
critical region of the nozzle decreases considerably 
with increased rate of injection. In the subcritical 
region, the decrease is much less pronounced. 

2. Wall injection moves the hydrodynamic throat down- 
stream from the physical throat. 

3. The calculated injected boundary layer thickness, 
based on a uniform plug-flow model, is approximately 
half as thick as the visualized layer. This is 
probably because in the calculations, an assumption 
was made that no mixing occurs between the injected 
and main stream fluids. 

4. For fixed stagnation conditions, the total flow rate 
through the nozzle is not altered significantly by 
injection. The mainstream inlet flow rate is dimin- 
ished by the amount of the injection upstream of the 
throat . 

5. Experimental results for both incompressible and com- 
pressible flow indicate that for blowing rates as 
large as 20$ of the ‘inlet mass flow rate, no serious 
disturbances occur in the flow field. The Froude and 
Mach number distributions are modified from the no- 
injection case, but the flow field remains smooth 
through the transcritical region of the nozzle. 

6. Schlieren pictures for the compressible nozzle flow 
indicate that it is difficult to obtain a thick blown 
layer at the nozzle throat unless significant injec- 
tion takes place in the upstream, subsonic region. 
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NOMENCLATURE 


A = flow area 
B s width of flow 

Cp = specific heat at constant pressure 

F = Froude number, V//gTT 

H = height of fluid 

h = enthalpy 

M = Mach number 

m = mass flow 

P = static pressure 

R - gas constant 

T s temperature 

V = flow velocity 

v = injection velocity 
w 

y = rat^o of specific heats of gas 
6 = blown layer thickness 
p = density 
0 = angle 

X = blowing parameter 


Subscripts 


inj = injected 

inlet = at inlet to nozzle (before injection) 
o = stagnation conditions 
t = nozzle throat 
w = at the wall 
« ~ at main stream conditions 
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FIGURE 2 MACH NUMBER AS A FUNCTION OF POSITION AND STAGNATION 
HEIGHT - NO INJECTION. 
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FIGURE 3 


SCHEMATIC OF NOZZLE SHOWING SYMBOLS AND CONTROL VOLUME. 
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FIGURE 4 EFFECT OF INJECTION ON POSITION OF HYDRODYNAMIC THROAT. 




FIGURE 5 FULL SCALE DRAWING OF NOZZLE SHOWING PQ:ROUS WALL SEGMENTS 
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FIGURE 12 EFFECT OF INJECTION? ON FROUDE NUMBER DISTRIBUTIONS . 
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FIGURE 13 PREDICTED FROUDE NUMBER AND HEIGHT DISTRIBUTIONS WITH 
INJECTION IN ALL THREE SECTIONS (v . = v TT = v IT } 
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FIGURE 14 PREDICTED FROUDE NUMBER DISTRIBUTIONS WITH 
INJECTION IN SECTIONS I AND II (v . = v 
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FIGURE 16 FULL SCALE DRAWING OF NOZZLE USED IN WIND TUNNEL. 



FIGURE 17 SCHLIEREN PHOTOGRAPHS OF AIR FLOW WITHOUT INJECTION. 



FIGURE 18 COMPOSITE PHOTOGRAPHS OF AIR FLOW IN WIND TUNNEL 
UPPER HALF - WITHOUT INJECTION 
LOWER HALF - WITH INJECTION OF AIR INTO AIR 
(WALL DENOTED BY LIGHT LINE) . 





FIGURE 19a AIR FLOW WITH LOW INJECTION RATES OF FREON IN 
ALL THREE SECTIONS. 



FIGURE 19b AIR FLOW WITH HIGH INJECTION RATES OF FREON 
IN ALL THREE SECTIONS. 




















IN FIRST SECTION ONLY. 




FIGURE 21 a AIR FLOW WITH HIGH INJECTION RATE OF 
FREON IN FIRST SECTION ONLY. 



FIGURE 21b AIR FLOW WITH HIGH INJECTION RATE OF 
FRFON IN SECOND SECTION ONLY. 




FIGURE 21c 


AIR FLOW WITH HIGH INJECTION RATES OF FREON 
IN FIRST AND SECOND SECTIONS. 








